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Abstract 

\o . 

. A model-independent analysis of decays of the form C — > £ + l~A (I = e, /i) is pre- 

sented, including the possibility that this three-body decay is preceded by an additional 
decay step D — > jC. Here A, C and D are heavy new-physics particles and j stands for 
a quark jet. It is assumed that A escapes direct detection in a collider experiment, so 
that one cannot kinematically reconstruct the momenta of the new particles. Instead, 
information about their properties can be obtained from invariant-mass distributions 
of the visible decay products, i. e. the di-lepton (U) and jet-lepton (j£) invariant-mass 
distributions. All possible spin configurations and renormalizable couplings of the new 
particles are considered, and explicit expressions for the invariant-mass distributions 
are derived, in a formulation that separates the coupling parameters from the spin and 
kinematic information. In a numerical analysis, it is shown how these properties can 
be determined independently from a fit to the mee and rrijt distributions. 



1 Introduction 



A large range of models have been proposed that predict new particles within the reach of 
the Large Hadron Collider (LHC). Since there is currently very little evidence for favoring 
one model over the others, it will be essential to analyze a potential new-physics signal in the 
LHC data in a model-independent approach, by independently determining the properties of 
each of the produced particles. Recently, this idea has gained increased interest, and several 
groups have worked on constructing such model-independent setups for a number of different 
observable signatures, see e.g. Refs. [IrH]. A particularly challenging scenario are processes 
that result in the production of new weakly interacting massive particles (WIMPs), which 
are invisible to the detector. WIMPs are predicted in many models as hypothetical dark 
matter candidates. In these models, the stability of the WIMP is a consequence of some 
(discrete) symmetry, under which it is charged. As a result, it can be produced only in 
pairs at colliders, leading to challenging events with at least two invisible objects. At hadron 
colliders like the LHC there are not enough kinematical constraints in events of this type for 
the direct reconstruction of the momenta of all particles involved. 

One approach to this problem is motivated by the fact that models predict additional new 
particles, which can decay into the stable WIMP. In this case, one can have cascade decay 
chains, which go through multiple decay steps before ending with the stable WIMP, so that 
one can construct invariant-mass distributions of the visible decay products of this cascade. 
The kinematic endpoints of these distributions yield information about the masses [5] of the 
new heavy particles, while the shape is sensitive to their spins [UEIEJE]- Refs. [HE] have 
analyzed decay chains built up from a sequence of two-body decays in a model-independent 
way, by considering arbitrary spin assignments [1] and also using general parametrizations 
for the coupling for the new particles [5]. 

However, for scenarios with relatively small splittings in the mass spectrum of the new- 
physics particles, it can often happen that the last decay step is a three-body decay mediated 
by a heavier off-shell particle, see right-hand side of Fig. [U In Ref. [8j, three-body decays have 
been analyzed in order to distinguish gluinos, the supersymmetric partners of gluons, from a 
Kaluza-Klein (KK) gluons in universal extra dimensions (UED) . A model-independent study 
of three-body decays has been presented in Ref. [9], but only in the limit of an asymptoti- 
cally large mass of the intermediate off-shell particles. In typical supersymmetry and UED 
scenarios, however, this limit is often not a good approximation. 

In this work, three-body decays of the form C — > £ + £~A will be analyzed in a model- 
independent setup without assumptions about the values of the masses of the new-physics 
particles. Here C is a massive new particle that decays into the WIMP A and two SM leptons 
£ ± = e 1 * 1 , u ± through the off-shell exchange of a third new particle B or the SM Z-boson, 
see Fig. [JJ. The spins of A, B, and C, their coupling parameters, and the mass m B of the 
particle B will be kept as free quantities that have to be extracted from the experimental 
data. We only impose the constraint tub > mc, or mz > mc — rriA, to ensure that we 

*In general, besides the Z-boson, a bosonic new-physics particle (e.g. a Z' or a Higgs boson) may also 
appear in the decay topology II. However, the branching of such a particle into leptons is strongly constrained 
by data on four-fermion contact interactions [10] . and thus its contribution will be neglected here. 
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Figure 1: Right: Three-body decays involving an off-shell new-physics particle B (topol- 
ogy I) or an off-shell Z boson (topology II). Left: The three body decay could occur as the 
last step of a longer decay chain. 

have an actual three-body decay. Without these constraints the three-body decay would 
decompose into two two-body decays, which is a scenario that has been discussed in detail 
in the literature cited above. 

Furthermore, we also consider the case that this three-body decay is the second step of a 

cascade decay of the form D — > q C — )■ q £ + £~A, where q refers to a SM quark (antiquark), 
see Fig. CD Such a decay chain would lead to two independent observable invariant-mass 
distributions, a di-lepton (£ + £~) invariant-mass distribution, and a jet-lepton (J^) invariant- 
mass distribution, where the jet emerges from the fragmentation of the quark or antiquark. 

For both of these cases, we investigate the simultaneous determination of the spins and 
couplings of the new particles A, B, C and D from the shapes of these distributions. The 
determination of the masses from kinematic endpoints has been discussed elsewhere [5], and 
here we will simply assume that the masses of the particles A, C and D are already known. 
On the other hand, the mass tub of the off-shell intermediate particle B can not be extracted 
from the kinematic endpoints, and we will study if instead it can be constrained from the 
shapes of the distributions. 

Our analysis closely follows the conventions of Ref. [2]. After introducing the relevant 
spin and coupling representations in section [2], the calculation of the 11 and j£ invariant-mass 
distributions is described in section[3j In sectionHJwe present a procedure for determining the 
spins and couplings of the new particles, as well as the mass of the intermediate particle B, 
by fitting the theoretically calculated functions to the experimentally observed distributions. 
The method is illustrated by applying it in two numerical examples. Finally, our main 
conclusions are summarized in section 
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2 Setup 



The three-body decay of a heavy new particle C into two opposite-sign same-flavor leptons 
and a second new particle A, 

C^ti'A, (£ = e,fi), (1) 

is mediated either by an off-shell heavy new particle B (with mg > mp > m^) or a SM 
Z-boson (with m c — rriA > rn z ). We also consider the possibility that eq. (pQ) occurs as the 
last step of a longer decay chain, 

D^qC 

U tCA. (2) 

Here D is a QCD triplet, while B and A/C are electrically charged and neutral QCD singlets, 
respectively. For the purpose of this work, it is assumed that A and C are self-conjugate 
(i. e. they are their own antiparticles)Q. Furthermore, it is assumed that A, B, C , and D are 
charged under some symmetry which ensures that A is stable and escapes from the detector 
without leaving a signal. 

In general, it is difficult to experimentally determine the overall strength of the couplings 
in the decay chain since the width of weakly decaying particles is typically small compared 
to the experimental resolution. Consequently, only the shape of the observable invariant- 
mass distributions will be considered here, similar to earlier studies on spin determination 
[TU2jl6f8jl9] • All expressions for these distributions presented in the following sections therefore 
include an arbitrary, but constant, normalization factor. 

Table [1] lists all possible spin assignments for the particles A—D in any renormalizable 
theory with fields of spin (scalars), spin 1/2 (fermions) and/or spin 1 (vector bosons). Also 
shown are examples for realizations of these assignments in known models. 

The chirality of the fermion couplings depend on the details of the new physics and 
thus are a priori unknown. Following Ref . [2] , we introduce arbitrary left- and right-handed 
components. For scalar-fermion-fermion vertices, the interaction Lagrangians are defined as 



■-A 



ip B A(a L uj- + a R u + )tp e + h.c, (3) 



B ~<Cr , ip A B (a L oj- + a R u + )ip e + h.c, (4) 



■-B 



ip c B (b L U- + b R u) + ) ip£ + h.c, (5) 



C — ~<r _ rl> B C(b L u- + b R fjj + )1> e + h.c, (6) 
~" B 



^Some new physics models predict decay chains with non-self-conjugate neutral heavy particles, which 
lead to distinct phenomenological features [11], but this case will not be considered here. 
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Table 1 : Possible spin configurations of the heavy particles D, C, B, and A in the decay 
chain of Fig.\l](F=Fermion, S=Scalar, V=Vector). Also shown are examples for realizations 
of these assignments in the Minimal Supersymmetric Standard Model (MSSM) or in models 
with one or two universal extra dimension (UED). Here q, £, and x° denote squark, slepton, 
and neutralino, respectively, q^, B® and W^'m refer to the first-level KK-excitations 
of quark, lepton, U(l) gauge Geld, and SU(2) gauge Geld, respectively. B° H ^ and W^-m are 
scalars stemming from one of the extra components of the higher- dimensional gauge Gelds 
in UED. More details of these models can be found in Refs. fl2 |, [T3 y . 



D "^Z^n ip D C (c L u)_ + c R u + )ip g + h.c, (7) 

^ q 

I > ip c D (c L u_ + c R u + )ij q + h.c., (8) 

where u± = |(1 ± 75). For vector-fermion-fermion couplings, A must be replaced by Jf± 
in ([3]), etc. After normalizing the overall coupling strength to unity, each vertex can be 
parametrized by a single angle a, (3, or 7, 

a L = cosa, &l = cos/3, c L = cos7, 
a R = sin a, b R = sin/3, c R = sin 7. 

As will be shown later, the entire parameter space for the couplings can be covered by 
restricting the angles to the intervals a G [— tt/2, 7r/2], /3,7 G [0, 7r/2]. 

The form of the CAZ vertices is uniquely determined by Lorentz symmetry and CP 
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properties (since the Z-boson is CP-odd, while the self-conjugate A and C are C-even): 

iCd^AZ", (10) 

-C^AZ", (11) 

-CApZ^ (12) 

{C»A V - ApOdPZ? + cy&., (13) 

iJcj^AZ", (14) 

where again the coupling constants have been normalized to unity. 

In an experimental analysis, it is impossible to tell on an event-by-event basis whether a 
quark or an antiquark is emitted in the first stage of eq. (T2]), i.e. whether the cascade decay 
was initiated by a particle D or its antiparticle D. However, the observable j£ invariant-mass 
distribution may depend significantly on the fraction / of events stemming from D decays 
versus the fraction / of events stemming from D decays, with f + f = 1. 

As pointed out in Ref. [2j, the ratio of / and / is very difficult to determine without 
model assumption and thus should be treated as a free parameter. The j£ distribution 
depends on / and / only through the combinations /|cl| 2 + /|cr| 2 = / cos 2 7 + / sin 2 7 and 
/|cr| 2 + J\cl\ 2 = f sin 2 7 + / cos 2 7. It is therefore convenient to introduce the parameter 7, 
defined by [2] 

cos 2 7 = / cos 2 7 + / sin 2 7, (15) 
sin 2 7 = / sin 2 7 + /cos 2 7. (16) 

From the analysis of the j£ invariant-mass distribution one can only obtain a constraint on 
7, but not on 7 and / independently. 




3 Invariant- mass distributions 

As pointed out above, it is difficult to discriminate experimentally between the decay chain in 
Fig. [TJ with a quark emitted in the first stage, and its charge-conjugated version with an anti- 
quark emitted in the first stage, since both quark and antiquark fragment into jets. Therefore 
the only relevant observable invariant-mass distributions are the mu (lepton-lepton) distri- 
bution and the rrijt (jet-lepton) distribution. 

There is no distinction between the two leptons in the three-body decay, in contrast to 
the situation when B can be produced on-shell (i.,e. for tub < rnc) in which case one can 
define a "near" and a "far" lepton [U[2l[5HZ]. 
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Explicit expressions for the m« and rriji distributions are obtained by computing the 
squared matrix elements for the different spin configurations 5=1—11 in Tab. [T] and inte- 
grating over the remaining phase space variables. A convenient choice for the phase space 
integration is given by 



1 dY 



T dm 2 u 



.max 

% al- 



N u / ' dm 2 M - \M 3 \\ (17) 



min.max 1 r 2 , 2 2-^x1/2/2 2 2 M 

m M 2 = 2 l m A + m c ~ m u T * {m A , m c , m u )\, 



r dm\ l+ 



r m cl 1 - m2 qe +/( m D-' m c)] r 27T 
= N qe dm 2 M - / dxj> (18) 



"(m 2 M _-m\)(m%-m 2 M _)/m? M _ ^ 

X / dm 2 u — 2 2— \M vl 



m C - m Ai~ 



4| ) 



where X(a, b, c) = a 2 + b 2 + c 2 — 2(ab + ac + be). In these equations, denote the 

matrix elements for the 3-body or (3+l)-body decay processes, respectively, while m^i- is 
the invariant mass of particle A and one of the leptons, and <fi is the angle between the plane 
spanned by the lepton-lepton system and the quark in the reference frame of C. The charge 
of the lepton in thai- an d m q t+ has been specified for definiteness, but one can equally well 
choose the variables rriAe+ and m q e-. Nu and N q i are unspecified normalization constants. 

The observable jet-lepton distribution dT/dm 2 e is obtained from dT/dm 2 e by replacing 7 
with 7, see eqs. (TT5 ]) . (fI5]) . 

As mentioned above, the endpoints of the invariant-mass distributions can be used to 
obtain information about the masses uia, rnc and mo of the particles that are produced 
on-shell in the cascade, while the shapes of the distributions depend on the couplings and 
spins of the particles A-D. Focusing on the latter, it is convenient to define the distributions 
dT /drhu and dT/drhj£ in terms of unit-normalized invariant masses 

™u = = m c - m A , (19) 

m je = ^ {mjn 2 = ~^{m 2 D - m 2 c )(m 2 c - m\). (20) 

For the spin configurations 5=1-6, the dependence on the coupling parameters a, /3, 7 can 
be cast into the form 

— = (cos 2 a sin 2 (3 + sin 2 a cos 2 0) f[ U \rh 2 u - : m 2 ,, m 2 B , m 2 c ) (21) 

r drriu 

+(cos 2 a cos 2 (5 + sin 2 a sin 2 0) f% (m 2 e ; m 2 Al m 2 B) m 2 c ) 
+(cosa sina cos/3 sin/3) ff i \rh 2 u \ m 2 A ,m 2 B ,m 2 c ), 
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r dm 



1 , 2 ■ 2 a 2~i-2 2 o • 2 ~\ Ajl) i ~ 2 2 2 2 2 \ / n n\ 

= (cos a sin p cos 7 + sin a cos p sin 7) /-f (m^; m A ,m B ,m c ,m D ) {22) 

1/ 2 ■ 2 o • 2 ~ 1 • 2 2 o 2~\ Ait) 1 ~ 2 2 2 2 2 \ 

+(cos asm p sin 7 + sm a cos p cos 7) j* ) (m je ; m A ,m B ,m c ,m D ) 

1 I 2 2 n 2~,-2 ■ 2 o • 2 ~\ Aj£)( " 2 2 2 2 2 \ 

+(cos a cos p cos 7 + sm a sm p sm 7) /g (mj t ; m A ,m B ,m c ,m D ) 

1/ 2 2 o • 2 - 1 • 2 • 2 o 2 ~\ Ajt) t " 2 2 2 2 2 \ 

+(cos ol cos p sm 7 + sm asm p cos 7) (mj t ] m A ,m B ,m ,m D ) 
+(cosa sin a cos/3 sin/3) f^ e \m^ e ] m\,m Bl m 2 Cl m 2 D ), 

where the functions fl and f± are independent of the coupling parameters, but they 
contain the entire kinematical and spin information, including the dependence on the par- 
ticle masses. Note that ff ] and receive contributions only from the interference term 
between the t- and ^-channel diagrams in the upper part of Fig. [Q see also Ref. [S]. 

From eqs. (|2ip . (|22|) one can see that without loss of generality the coupling parameters 
can be restricted to the intervals a G [— tt/2, 7r/2], /3,7 G [0,7r/2], as already mentioned in 
the previous section. 

For 5=7-11, the CAZ coupling is uniquely fixed up to an overall coupling constant, so 
that there is only one term for the lepton-lepton invariant-mass distribution. However, there 
are two possible terms for the jet-lepton invariant-mass distribution: 

T^ = f W (^rn%m%,ml), (23) 

1 di\- „ ( *j f \ /O Q Q Q Q \ n( *} P\ /O O O Q Q \ /\ 

P dft~ = fs \ m j6 m A ,m z ,m c ,m D ) + cos 2j fX "{m je ; m A ,m z ,m c ,m D ), (24) 

The lepton-lepton distribution dT/drhu can be expressed in terms of compact analytical 
formulae. On the other hand, the analytical results for dT/dm q e are very lengthy, so that 
instead we chose to perform the last integration step (over m 2 M „) numerically. 

Explicit expressions for the functions are available for free download (see appendix). 
Figs. [2H31 depict the distribution functions for a sample mass spectrum. In the figures, the 
overall normalization constants have been fixed by requiring that f\ , f\ , f^ u \ and /g"' 7 '* 
are unit-normalized. The right column of Fig. [2] also illustrates how the distributions vary 
with the mass of the off-shell intermediate particle B, for the example of the spin 
configuration 5=1. 



4 Analysis method 

In this section we will discuss the determination of the spins and couplings parameters of 
the new particles, as well as the mass of the off-shell particle B, by fitting the theoretically 
calculated distributions to experimental data. The general procedure will be outlined in the 
next subsection, while its application will be demonstrated in subsection I4.2l for two concrete 
numerical examples. 
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Figure 2: Left: Distribution functions (i = 1,...,3) for the spin configurations S=l-6, 
for tub = 200 GeV. Right: Dependence of f± (i = 1,...,3) on the mass ms of the intermedi- 
ate particle for the case S=l. The other mass parameters have been chosen as mc = 184 GeV 
and rriA = 98 GeV. In these plots the overall normalization has been fixed by normalizing 
fi to unity. 
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Figure 3: Distribution functions (i = 1,...,5) for the spin configurations S—l-6. The 
mass parameters have been chosen as vtld = 565 GeV, mc =184 GeV, mg = 200 GeV and 
rriA = 98 GeV. In these plots the overall normalization has been fixed by normalizing 
to unity. 
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Figure 4: Distribution functions and /g A for the spin configurations S=7-ll. The 
mass parameters have been chosen as vhd = 565 GeV, mc = 184 GeV and tua = 98 GeV. In 
these plots the overall normalization has been fixed by normalizing and f^ e ^ to unity. 



4.1 Conceptual procedure 

The analysis is based on a binned x 2 fit for the ii and ji distributions. In this fit, the binned 
histogram for the data is compared with theoretical histograms obtained by numerically 
integrating the functions f± and f i , defined in the previous section, over the interval 
covered by each bin. In the fit, the coupling parameters a, (3, 7 and the mass m B are kept 
as free parameters. Varying over these parameters and the spin configuration S, the best-fit 
result is found as the set of numbers {S, a, (3, 7, mg} that minimizes the \ 2 value. 

During the fit procedure, for every given choice of the parameters {S, a, (3, 7, m B }, the 
theoretical histograms for the ii and ji distributions are normalized such that the total 
number of events in the theoretical histogram agrees with the number of events in the data 
histogram. In practice, this normalization is most easily carried out numerically. 

In general, it may happen that there is not a unique solution for the minimum x 2 value, 
but instead several degenerate best-fit points are obtained. In such a situation, the coupling 
parameters a, (3, 7 and/or the spin assignment S cannot be determined uniquely from the 
observable distributions of the decays ([I]),© alone. 



10 



4.2 Numerical examples 

To illustrate the fitting procedure, its application is demonstrated by performing a fit to 
mock-up data histograms. This section is based on the parton-level description of the de- 
cay processes ([I]),© as described in the previous sections, thus neglecting issues such as 
backgrounds, jet combinatorics and energy smearing, which are relevant in a realistic exper- 
imental setup. However, earlier studies [SJIH] have shown that, for mass parameters similar 
to the ones chosen here, it is possible to obtain a clean, almost background-free sample of 
signal events with relatively simple selection cuts. 

Let us consider two sample choices for the hypothetical data: 

"Data" A: 5 = 1, a = 0, = vr/2, 7 = 0, 

m D = 565 GeV, m c = 184 GeV, m B = 200 GeV, m A = 98 GeV 

(corresponding to the MSSM decay chain qi — > X2 ~~ ^ ~ > 

"Data" B: 5 = 11, 7 = 0, 

m D = 565 GeV, m c = 184 GeV, m A = 98 GeV 

(corresponding to the MSSM decay chain q L — > X2 ~ * Xi)- 

For each case, we have computed "data" histograms with 10 bins each for the mm and the 
rhj£ distributions, corresponding to a total of 1000 events. Then we have performed a \ 2 fit 
of the theoretical distribution functions to these fake "data" histogram for each of the spin 
configurations 5=1-11, searching for the minimum x 2 value as a function of the parameters 
a, f3, 7, and mj§. 

The results are shown in Tables |2] and EJ From Tab. [2] one can see that when only 
information about the mu distribution is available, it is difficult to distinguish the "data" A 
(based on the spin configuration S=l) from the spin configurations 5=2-6. The underlying 
reason is that for each of these spin configurations there are three unknown continuous 
parameters, a, (3 and ms, which can be adjusted so as to mimic the data distribution. 

On the other hand, the spin configurations 5=7-11 can be distinguished from "data" 
A with high significance, using only the m« distribution. This is a consequence of the fact 
that there are no free parameters to adjust in dT/drhu for 5=7-11, and that these spin 
configurations correspond to a different diagram topology (Topology II in Fig. [1] instead of 
topology I). 

If both the mu and rhj£ distributions are included in the fit, all possible spin configurations 
can be discriminated with at least six standard deviations, for the given number of 1000 
events. 

For the second example, it is evident from Tab. |3] that "data" B can be distinguished 
from all other spin configurations 5=1-10 by just using the mu distribution. In fact, for all 
combinations except 5=8 and 5=9 the significance for this discrimination is very high and 
is not improved substantially by including the rhji distribution in the fit. Also note that the 
best-fit results for 5=1-6 are obtained for very large values of mg, since increasing values 

••■For the spin configurations 5=7— 11, the non-zero Z-boson width has been included although its numerical 
impact is not very important for the masses chosen here. 
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a) "Data" A, using only rhu distribution: 



s 


min x 2 


best-nt naramptprs 
a (3 iris [GeV] 


1 [SFSF] 


0.00 


0.00 1.57 200.0 


2 [FSFS] 


0.00 


-1.22 1.05 209 


3 [FSFV] 


0.00 


+1.14 0.43 197.7 


4 [FVFS] 


0.27 


-1.34 0.23 216 


5 [FVFV] 


0.05 


-0.38 0.38 197 


6 [SFVF] 


0.05 


-0.65 0.92 191.3 



b) "Data" A, using both rhu and rhji distributions 



s 


min x 2 


bpst-fit paramptprs 
a (3 7 m B [GeV] 


1 [SFSF] 





0.00 1.57 0.00 200.0 


2 [FSFS] 


150 


-0.08 0.07 1.57 754 


3 [FSFV] 


87 


±1.57 1.57 0.29 210 


4 [FVFS] 


48 


±1.19 0.00 1.57 220 


5 [FVFV] 


46 


-0.93 0.25 1.57 224 


6 [SFVF] 


37 


-0.50 0.53 1.57 197.4 



5 


min x 2 


7 [FSS] 


140 


8 [FSV] 


3100 


9 [FVS] 


4200 


10 [FVV] 


290 


11 [SFF] 


3700 



5 


min x 2 


best-fit 

7 


7 [FSS] 


200 


? 


8 [FSV] 


3100 


? 


9 [FVS] 


4300 


0.39 


10 [FVV] 


330 


1.57 


11 [SFF] 


3700 


1.08 



Table 2: Results for fitting all spin configurations S—l-11 to (a) the rhu distribution only, 
and (b) the rhu and rhji distributions together, using scenario "data" A for the mock-up 
data histograms. Shown are the minimum x 2 (rounded to two significant digits) for each 
spin conhguration, as well as the parameter values for which this minimal value is attained. 
"?" indicates that the x 2 value is independent of that parameter. The numbers correspond 
to 1000 events. 



of rriB shift the rhu distribution toward larger values of rhu, see Fig. [2] (right), leading to 
better agreement with the reference case 5=11, see FigHl 

In addition to the spin determination, the couplings of the new particles and the mass 
of the off-shell B particle can in principle be extracted from the fit to the invariant-mass 
distributions. This is shown in Fig. for the example of "data" A. The panels (a) and (b) in 
the figure depict the constraints on a, (3 and m B obtained from fitting the rhu distribution 
alone, assuming that 5=1 is the correct spin configuration. If a fit to both the rhu and rhji 
distributions is performed, one obtains the results in panels (c) and (d). As evident from the 
plots, the inclusion of the rhji distribution does not only lead to a constraint on 7 (which 
cannot be obtained from dT/drhu), but also to improved bounds on a and (3. 

However, the fit results for the coupling parameters always have a two-fold degeneracy, 
since the invariant-mass distributions, eqs. (!2T]) - (!24l) . are invariant under the transformation 
{a, (3, 7} -+ {signa (f - \a\), | - p, | - 7}. 
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a) "Data" B, using only rhu. distribution: 



s 


min x 2 


best-nt naramptprs 
a (3 rriB [GeV] 


1 [SFSF] 


1200 


+0.79 0.79 oo 


2 [FSFS] 


670 


? ? oo 


3 [FSFV] 


2200 


? ? oo 


4 [FVFS] 


1100 


? ? oo 


5 [FVFV] 


720 


a = (3 = 7 oo 


6 [SFVF] 


740 


±1.57 0.00 oo 



b) "Data" B, using both rhu and rhj£ distributions 



s 


min x 2 


best-fit parameters 
a (3 7 m B [GeV] 


1 [SFSF] 


1200 


±0.78 0.77 0.00 oo 


2 [FSFS] 


690 


? ? ? oo 


3 [FSFV] 


2300 


? ? ? oo 


4 [FVFS] 


1100 


±1.25 0.43 1.32 oo 


5 [FVFV] 


750 


±0.46 0.46 1.57 oo 


6 [SFVF] 


760 


±1.57 0.00 ? oo 



S 


min x 2 


7 [FSS] 


1600 


8 [FSV] 


16 


9 [FVS] 


8.7 


10 [FVV] 


1100 


11 [SFF] 


0.00 



s 


min x 2 


best-fit 

7 


7 [FSS] 


1600 


? 


8 [FSV] 


25 


? 


9 [FVS] 


59 


0.00 


10 [FVV] 


1100 


0.00 


11 [SFF] 


0.00 


0.00 



Table 3: Same as Fig. \2\ but using "data" B for the mock-up data histograms. 



5 Summary 

In this paper, a general analysis of three-body decays of the form C — > £ + £~A, leading to 
a pair of opposite-sign leptons and one invisible particle A, has been presented. This decay 
process can occur in many proposed new-physics models, either from direct production of 

the particle C at the LHC, or from a cascade decay of the type D — > <1 C — )■ <? £ + £~A, both 
of which have been studied here. 

No assumptions about the masses, spins and couplings of the participating new-physics 
particles have been made, including the off-shell particle B mediating the three-body decay. 
Instead, all possible spin configurations and coupling form factors have been considered. Ex- 
perimentally, the masses, spins and coupling parameters may be determined from measuring 
the invariant-mass distributions of the visible decay products. 

In the present case, there are two independent distributions, one with respect to the 
di-lepton (£ + £~) invariant mass, and the other with respect to the jet-lepton (j£ ± ) invariant 
mass. Results for both have been obtained in terms of relatively compact analytical functions 
or one-dimensional integral representations. 

In two concrete numerical examples, it has been tested how well the properties of the 
new-physics particles A, B, C and D can be determined from these two invariant-mass 
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tt/6 tt/3 nil 190 195 200 205 210 x/6 n/3 n/2 190 195 200 205 210 



Figure 5: Determination of the parameters a, /3, 7, and vtib using only the rhgi distribution 
(a,b), and using both the m« and rhjg distributions (c,d). The dark/light bands in the figures 
correspond to the 68%/95% confidence-level regions. The plots correspond to a sample of 
1000 events for the scenario "Data" A. 

distributions. It turns out that the di-lepton invariant-mass distributions alone is sometimes 
not sufficient to uniquely determine the spins and coupling parameters. However, if the 
longer two-step cascade decay chain is observed, and one can measure both the £ + £~ and jt^ 
invariant-mass distributions, it is possible to unambiguously discriminate between all possible 
spin configuration with high significance. Furthermore, one can independently constrain all 
coupling parameters and the mass of the off-shell mediator B, up to an intrinsic two-fold 
ambiguity. 

The results presented here are based on a parton-level analysis. In a realistic experi- 
mental environment, the significance for the model discrimination and the precision for the 
parameter determination may be diluted by jet energy smearing and combinatorics, but the 
essential features and main conclusions are not affected substantially by these effects. 
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Appendix: Formulae for invariant- mass distributions 



Explicit expressions for the functions f\ and f> are avaiable in Mathematica format at 
http://www.pitt.edu/~afreitas/dec3.tgz. Note that the expressions in this file are not 
normalized, since in practice the normalization is best carried out numerically as described 
in section 14.11 The results for ff^ are given as analytical formulae, while are presented 
in terms of one-dimensional integral representations of the form 

fM= dm\,.Ff\ (25) 

Jm 2 A 
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